Local modulations of the spin-fluctuation mediated pairing interaction by impurities 

in d-wave superconductors 
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We present a self-consistent real space formulation of spin-fluctuation mediated d-wave pairing. 
By calculating all relevant inhomogeneous spin and charge susceptibilities in real space within the 
random phase approximation (RPA) , we obtain the effective pairing interaction and study its spatial 
dependence near both local potential and hopping impurities. A remarkably large enhancement of 
the pairing interaction may be obtained near the impurity site. We discuss the relevance of our 
result to inhomogeneities observed by scanning tunneling spectroscopy on the surface of cuprate 
superconductors. 

PACS numbers: 74.20.-z, 74.25.-q, 74.62. En, 74.72.-h 

I. INTRODUCTION 

Scanning tunneling spectroscopy (STS) experiments 
on surfaces of several high-T c materials* have discov- 
ered that the electronic structure of at least some of the 
cuprate superconductors is strongly inhomogeneous in 
real space. Modulations observed include checkerboard 
local density of states (LDOS) patterns which peak at 
biases near the local gap energy, together with strong 
gap modulations on nanometer length scales^— While 
the latter modulations appear to be random in charac- 
ter and driven by disorder, the STS conductance maps 
contain important correlations of gaps with other phys- 
ical observables. For example, the size of the local gap 
was found to be positively correlated with atomic scale 
defects, thought to be interstitial oxygen dopants^ (see 
however Ref. |8|). This observation inspired a reexamina- 
tion of the effect of an impurity in these systems, since 
it had been anticipated that an oxygen, which donates 
two holes to the CuC>2 plane, would locally overdope the 
system and lead to a smaller gap in the neighborhood of 
such a dopant impurity. The positive correlation between 
the putative oxygens and the observed spectral gap led 
Nunner et alM-^ to propose that the dopants might be 
enhancing the pair interaction on the atomic scale. This 
might occur if the local electronic structure was modified 
significantly and seems to be consistent with a number 
of other recent STS results 

Nunner et ai£ observed on phenomenological grounds 
that this simple assumption of dopant atoms modulat- 
ing the pair interaction could explain a large number 
of observed correlations; for example, it correctly repro- 
duces the anti-correlation of coherence peak height and 
peak energy, the correlation of impurity position with 
gap size, and the detailed frequency dependence of the 
0:LDOS(w) correlation. In almost all treatments of dis- 



order in superconductors, impurities are assumed to sim- 
ply scatter electrons via a screened Coulomb potential, 
rather than modify the pair interaction locally. The well- 
known suppression of the gap due to pairbreaking is of 
course observed in such theories if they treat the gap 
self-consistently; this effect is distinct, however, from a 
modulation of the underlying interaction, which may be 
expected to influence the gap in a potentially new way. 
Exceptions to the standard approach in the literature 
on conventional superconductors include Refs. [l6l where 
the notion of a local pair interaction was considered in 
some specialized contexts; from these works it is clear 
that these effects will be largest in systems where the 
coherence length is small, like the cuprates. 

The goal of such theories has been to use the STS data 
to potentially deduce which changes in electronic struc- 
ture influence the pairing interaction locally, and thereby 
learn about the origin of pairing itself in these systems. 
For such a program to be successful, one needs an under- 
standing of the location of the dopants and the changes in 
electronic structure which they cause, as well as a micro- 
scopic model capable of including these changes. Thus far 
some work has been performed on the first question. For 
example, He et aliU performed density functional the- 
ory (DFT) supercell calculations to find the position of 
lowest energy for an O dopant in a E^S^CaC^Osunit 
cell, and noted that the interstitial atom caused a rigid 
tilting of the Cu04 half-octahedron, as well as giving rise 
to an additional weakly hybridized state at -leV which 
coupled well to the STM tip. 

In addition, considerable model-based work has been 
done on the second question, that of determining how an 
impurity might actually affect the pairing interaction. In 
the strong coupling limit, superconductivity is sometimes 
treated within slave boson mean field theories of the t — J 
model ) 18 ' 19 whereby the pairing interaction is provided 
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by the super-exchange interaction Jy itself. Assuming 
some random distribution of exchange couplings 
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sults were obtained for the gap maps and other observ- 
ables which were similar to those found by Nunner et 
al£ The study of the microscopic origin of pairing inter- 
action modulations was begun using the strong coupling 
approach by Maska et al.^ who assumed that the pres- 
ence of the dopant atoms induces a position-dependent 
shift of the atomic levels in the CuC>2 plane, and then 
calculated the locally modified super-exchange interac- 
tion J between copper ions in the presence of an impurity 
from a perturbation expansion of the one-band Hubbard 
Hamiltonian up to second-order. Within this framework, 
the diagonal disorder in the plane always leads to an en- 
hancement of J; accordingly, with the assumption that 
pairing is due to super-exchange, the pairing supercon- 
ducting gap always increases in the vicinity of the dopant 
atoms, as assumed phenomenologically in Ref. [^. How- 
ever, Johnston et al£^ argued based on cluster model 
calculations for the three-band Hubbard Hamiltonian^ 
which account explicitly for the Cu-0 hopping processes, 
that electrostatic modifications due to the presence of 
the oxygen dopants locally suppress J. An enhancement 
of J was only possible if they assumed that electronic 
coupling to local phonon modes was strongly modified 
by the dopants. Foyevstova et al~ returned to the per- 
turbation theory approach in the three-band model, and 
showed explicitly that whether an impurity enhances or 
suppresses the local exchange in this somewhat more real- 
istic case depends explicitly on details of the Cu-0 model 
parameters. 

For completeness it should be mentioned that a re- 
cent calculation within the Gutzwillcr approximation to 
the t — J model including disorder also studied the lo- 
cal disorder-induced variations in the effective super- 
exchange couplings J.— This study, however, focused on 
the consequences of charge reorganization and the doping 
dependence of disorder-induced magnetism, and did not 
include the explicit J-enhancements caused by renormal- 
ized local Coulomb interactions.— 

Assuming, as in the strong coupling case, that the 
mean field interaction can be taken as a local pairing 
strength is intuitively reasonable, but the approxima- 
tion is even less well controlled than in the homogeneous 
case. It should be possible to study this problem from 
the weak coupling side as well, in an approach where 
the effective interaction in the particle-particle channel 
can be clearly defined, and which has the additional ad- 
vantage that it can be more easily coupled to ab ini- 
tio results which provide the needed structural and elec- 
tronic structural details for a given system. An "aver- 
age" version of the calculation envisaged was performed 
by Foyevstova et aZ.j2£ who calculated the total singlet 
spin-fluctuation interaction^ for a Bi2Sr2CaCu20s sys- 
tem with and without an O dopant; in this case the defect 
was indeed found to enhance the overall pairing. 

Nevertheless, one would like to understand the influ- 
ence of an impurity on pairing at the local level. In this 



work we continue this program from the microscopic side 
by formulating the full RPA particle-particle vertex in 
the singlet channel in real space for a general inhomo- 
geneous system. This interaction has a somewhat more 
complicated mathematical form compared to the popu- 
lar Berk-Schricffcr fluctuation exchange expression] 29 ! 30 
We then evaluate the real pairing interaction for a single 
impurity in a Hubbard model, and show that remarkably 
large enhancements of pairing strength can be obtained 
near a dopant impurity. The spin-fluctuation model pro- 
duces a characteristic shape of the pairing modulation 
near a strong impurity which is reflected in the LDOS 
pattern observable by STM. In our conclusions we dis- 
cuss the prospects for observing this behavior. 



II. MODEL AND METHOD 

In the initial step of the calculation of the effective 
real space pairing potential, we obtain the electronic den- 
sities calculated self-consistcntly in the normal state us- 
ing a mean-field approximation to the one-band Hubbard 
model 

H = ^ Ujc^Cja + (n ta ) - p)n lS 

J~] VimpSjn - r imp )n i(T . (1) 



Here, c\ a refers to creation of an electron with spin a 
at lattice site i, and rii a is the number operator of spin 
<7 particles at site i. The hoppings ti j are included to 
nearest t = 1, and next-nearest neighbor sites t' = —0.3. 
Note that this Hamiltonian also contains the impurity 
potential V; mp at a site placed at position r^. . A di- 
agonalization of Eq.(fT]) allows us to obtain the effective 
interaction V e s(i,j) which in real space can be written as 



V eB (i,j) = U+- s 
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The effective interaction @ is due to longitudinal and 
transverse spin-fluctuations using the approach of Berk 
and Schrieffer^ and the derivation is given in the Ap- 
pendix. The susceptibilities entering Eq.([2]) are real 
space matrices and given by 



f(E m cr) — f(E n a') 



En 



E„ 



II] 



(3) 



in terms of the eigenvectors u ma and eigenvalues E ma 
obtained in the diagonalization of Eq. (fTJ) . 

It is useful now to consider the partial Fourier trans- 
form of the susceptibility with respect to relative coor- 
dinate r = Vi — Tj, leaving explicit dependence on one 
spatial variable, 



x(q,r t ) = J2> 



.iq^ri-r,-) 



(4) 
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Note that Eq. <j2j) reduces to the usual form in the case 
of translational invariance, x( c L7 r i) = I n addition, 

it is worth remarking that the convention used here is 
slightly different from other works, e.g. Ref. HH where 
the mixed susceptibility is a function of the center of mass 
variable (r; + ry)/2. 

After the calculation of the effective spin-fluctuation 
mediated pairing, the densities (nj CT ) and superconduct- 
ing gap values Ajj are calculated sclf-consistently from 
Eq.([T]) with the addition of a standard BCS singlet pair- 
ing term, i.e. from the Hamiltonian 



so 



hJ 



where Aj i3 - = Y^^Jl (c^c^ — CjfCii). The factors of 1/2 
arises from the restriction to the singlet pairing channel, 
cf. the Appendix. 
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III. RESULTS 

A. Point-like nonmagnetic impurity 

In Fig. [TJa) we show the normal state electron den- 
sity, n, = (n^ + nil) as a function of site i around 
a single nonmagnetic point-like impurity. In contrast 
to the superconducting case, in the normal state there 
is in this model never any induced local magnetization 
due to the absence of the impurity bound state near 
the Fermi level.— ~— The effective pairing interaction be- 
tween the nearest neighbor sites is shown in Fig. [TJb), 
where we at each site show the average interaction to 
the four nearest neighbors: V sc (i) = j^2jVes(i,j), 
j 6 {i ± x, i ± y}. It is seen that the effective pair- 
ing interaction is suppressed at the impurity site, but 
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FIG. 1: (a) Local charge density m in the normal phase 
prior to calculation of the effective pairing interaction. For 
the results presented here the parameters are: U = 2.2, 
t' — —0.3 (all energies are given in units of t) and doping 
x — 0.15. A strong point-like nonmagnetic impurity 
(Fi m p = 10) is situated at site (xi mp , 2Am P ) = (13, 13). (b) 
Effective pairing interaction between nearest neighbors for 
the same system as in (a). For each site the average 
potential to the four nearest neighbors is plotted. 



FIG. 2: (a) Local charge density m in a system with a 
strong nonmagnetic impurity (Vi mp = 10) and a spatially 
homogeneous superconducting pairing potential 
V e ff = —0.43. Parameters: U = 2.2, t' — —0.3 and doping 
x — 0.15. (b) Same as (a) but in the case of the spatially 
inhomogeneous spin-fluctuation mediated pairing potential 
V e fi(i,j) from Fig. [TJb). (c,d) Local superconducting d-wave 
order parameter corresponding to the two cases in (a) and 
(b), respectively. (e,f) Local density of states at energy 
E = 0.78 (close to maximum gap value) corresponding to 
the two cases in (a) and (b), respectively. 



enhanced on a quartet of small regions close to the im- 
purity. The response of the system to a perturbation is 
dominated by the peaks in the homogeneous susceptibil- 
ity, although because the perturbation is space dependent 
all wave vectors are coupled. In the present band struc- 
ture, x RPA (q,^) = X°(q,w)/(1 - U X °(q,oj)) at u = is 
peaked at incommensurate wave vectors Q + Aq, where 
Q = (7r,7r) and Aq ~ (±0.5tt, 0), (0, ±0.5tt). The spa- 
tial modulation of the nearest neighbor pairing shown in 
Fig. [TJb) is determined by Aq, giving rise to the quartet 
of peaks at roughly four lattice spacings from the im- 
purity site. With the effective pairing interaction given 
in Eq.©, we can calculate self-consistently the densities 
(riia) and the superconducting d-wave order parameter 
Aij from Hgc, sec Eq.©. For all results shown in this 
paper we include the twelve closest sites j around each 
site i in V e fi(i,j). The resulting spatial structure of the 
density modulations and the superconducting gap varia- 
tions are shown in Fig. [2^b,d). The superconducting or- 
der parameter is averaged over the four nearest neighbor 
sites Aj = i[A. M+4 + A. M _ S - A^+g - A. M _ S ]. For com- 
parison, we also calculated the densities and the d-wave 
order parameter in a system containing the same impu- 
rity, but with a homogeneous pairing potential calculated 
from a clean system with the same parameters, giving a 
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FIG. 3: Maximum of the superconducting d-wave gap, 
A ma x; as a function of Coulomb interaction U for point-like 
nonmagnetic impurities and a homogeneous system. For all 
runs the system displays no antiferromagnetic order, since 
the critical U in this system is U C 2 > 2.4. The inset shows 
the profile of the d-wave gap as a function of site (x,yi mp ) 
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FIG. 4: Average of U X n(Q) around Q = (ir, n) for different 
values of U as a function of site Vi = (a:,j/i mp ). Black lines 
are C/(xtJ.(Q)) f° r a homogeneous system whereas colored 
lines are for systems with a nonmagnetic impurity at 
(Ximp, 2/imp) = (13, 13). 



nearest neighbor pairing interaction V^ c = —0.43. The 
result is shown in Fig. &[c) , where it is seen that no en- 
hancement of the d-wave gap is present, as opposed to 
the case of spin- fluctuation mediated pairing (Fig. HJd)), 
where the d-wave gap is locally enhanced and the spa- 
tial structure is similar to the structure of the effective 
pairing potential from Fig. HJb) . Note that the scales of 
the pairing interaction V e s and superconducting gap A 
modulations are similar. In a clean superconductor, the 
order parameter response to a local perturbation will be 
smeared out over a coherence length £oj which in a re- 
alistic system is larger than the atomic lattice constant 
a; here, however we have shown results for an artificially 
large gap, corresponding to a coherence length £o ~ a, in 
order to display the desired physical effects more clearly. 

The LDOS is suppressed at the impurity site both in 
the case of a constant effective interaction and in the 
case of spin-fluctuation mediated pairing, see Fig.[2{c,f). 
However, whereas the LDOS in the case of constant ef- 
fective interaction is roughly constant for all sites away 
from the impurity the picture is different when the pair- 
ing interaction is due to spin fluctuations. In the latter 
case, we find significant LDOS enhancements at energies 
close to the gap value at the same quartet of real space 
regions seen in Figs.QJb) and^d). 

From Fig. ^d) it is clear that certain sites close 
to the impurity site exhibit a local maximum of the 
superconducting gap. In Fig. [3] we show how this 
maximum gap value depend on both the strength of the 
Coulomb interaction U and the impurity potential Vi mp . 
For reference, we have plotted the gap in a homogeneous 
system which is also enhanced upon increased U. This 
effect arises simply from the singularity in the RPA 
susceptibility. Since we tune the system towards the 
singularity from below in the paramagnetic phase, it is 
the last term in the effective interaction of Eq.([2]) which 



is the most important. In the system containing a single 
impurity, the gap is strongly enhanced as was shown in 
Fig. |3l The local structure of the gap enhancement is 
similar for a strong and a weak impurity as seen from the 
inset in Fig. [3] which shows a profile of the d-wave gap 
through the impurity site. However, the enhancement is 
much more pronounced in the strong impurity case. 
The local enhancement effect can be understood from 
spatially-dependent variations in the bare transverse 
spin susceptibility. Due to the finite value of next-nearest 
neighbor interaction, t' — —0.3, the susceptibility dis- 
plays splitting in reciprocal space into incommensurate 
peaks and a broad maximum around the ordering vector 
Q = (tt, 7r). We Fourier transform the real space suscep- 
tibility to obtain xtl( ( li r i) an d map out the change in 
the magnitude of the susceptibility near Q as a function 
of site r,. To account for the broadening of the peak and 
to minimize finite-size effects, we average over a region 
containing Q and the eight closest g-values around Q 
to obtain (x°.(Qj r i))- The result of this procedure is 
given in Fig. [4] It is seen that due to local variations 
in the bare transverse spin susceptibility, the quartet of 
points with enhanced pairing in Figs.QJb) and[2Id,f) are 
characterized by being locally closer to the Stoner insta- 
bility Ux®i (Q) — ► 1 than any other sites in the system. 
This is the origin of the observed enhancement of the 
pairing potential and ultimately what causes the local 
gap enhancements. Specifically, charge density varia- 
tions due to an impurity act to tune the system closer 
to the Stoner instability locally which results in gap 
enhancements in the neighboring regions of the impurity. 
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FIG. 5: (a,b) Effective pairing interaction between nearest 
neighbors. For the results presented here the parameters 
are: U = 2.2, t' = —0.3, doping x — 0.15. A reduced 
hopping impurity of strength —20 % (a) and —100 % (b) is 
situated at site (xi m p,yi mp ) = (13,13). (c,d) Local 
superconducting ei-wave order parameter for the two 
impurity cases. (e,f) Average of E/x?4.(Q) around Q = (tt,ty) 
as a function of site r; = (rr,j/i mp ). 



B. Impurity with reduced hopping constants 

In this section we present the results of the local pairing 
modulations for another kind of impurity characterized 
by all hopping constants onto the impurity site being 
reduced by a certain percentage. In the extreme limit 
of 100% reduction, we obtain results similar to the case 
of a strong nonmagnetic potential scatterer as expected 
(compare e.g. Fig. [2jd) and Fig. [5jd)). Both magnitude 
and local structure of the superconducting d-wave gap 
is similar in these two cases. Furthermore, there is no 
induced local antiferromagnetic order. In the interme- 
diate regime, however, where the reduction is less than 
100%, we find a different local structure of the <i-wave 
gap enhancement. As an example, we show in Fig.[5ja,c) 
the effect of reducing the hopping elements by 20%. In 
this case weak local antiferromagnctism is induced (not 
shown), and the enhancement of the local pairing inter- 
action as well as the d-wave gap is maximal at the im- 
purity site. Again this gap structure can be understood 
from the local modulations in the transverse spin suscep- 
tibility, see Fig. [5{e) . Note that for the effective pairing 
interaction shown in Fig.(5Ja), the spatial structure of the 
effective interaction and the size of the amplitude mod- 



ulations are closely comparable to the phcnomenological 
values used in Ref. IH to model the STS gap variations in 
optimally doped Bi 2 Sr 2 CaCu208. 



IV. CONCLUSIONS 

The results of the local gap enhancements found in 
this paper gives rise to a number of additional questions 
for future studies. First, is the LDOS within the present 
formalism from more realistic many-impurity simulations 
consistent with the experimental STS data? As shown by 
Nunner et al£ the positive correlation between dopant 
atoms and large gap regions will certainly exist with the 
presented scenario, but the LDOS evolution with dop- 
ing, where pseudogap regions become dominant in the 
undcrdopcd regime cannot be captured within this weak- 
coupling approach. In addition, several unsolved LDOS 
issues remain regarding the origin of the so-called extinc- 
tion of impurity resonances near strong onsite impurity 
potentials in large-gap regions.— 

Second, what are the effects on a modified pairing 
near other kinds of impurities? We have recently applied 
the present formalism to also study the pairing inter- 
action and local gap structure near point-like magnetic 
impurities^ There it was found that the local structure 
of the gap enhancement is very similar compared to the 
case of point-like nonmagnetic impurities even though 
local antiferromagnetism is induced around a magnetic 
impurity site. 

Third, what is the feedback effect on the pairing inter- 
action in coexisting phases where e.g. spin- and charge 
density waves are present? The current formulation of 
the real space pairing mechanism does not hold in the 
ordered state. The simple expression that was derived 
in Eq.([2]) predicts a singularity in the effective pairing 
interaction as the system develops long range antiferro- 
magnetic order, as has been discussed previously. This 
effective vertex is, however, known to be suppressed by 
interactions with spin waves (Adler principl o 36 i 37 ), and 
our model is therefore only valid in the paramagnetic 
phase. 

Fourth, what are the thermodynamic consequences of 
the local enhancements of the gap found above? In par- 
ticular one might expect a non-trivial evolution of the 
residual resistivity Ap and T c suppression with increasing 
disorder concentrations. The slow rate of T c suppression 
relative to Ap in the case of strong scatterers such as Zn 
in cuprates has never been satisfactorily explained, but 
has been attributed to correlation effects 4£ This inter- 
esting question will be the topic of a future study. 

Finally, how robust is the present result of local en- 
hancement of pairing to the form of the bare interaction 
vertex? In the strong coupling limit, Foyevstova et al££- 
showed that the local pairing interaction (exchange) was 
increased only in the limit of pure Hubbard-type corre- 
lations in a single band^ 3 - We have established here the 
weak coupling analog of this result, but it is clear that 
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other types of interactions in the host, beyond the scope 
of this work, might lead to different results. 

Summarizing, we have studied how single nonmagnetic 
potential and hopping impurities modify the local pair- 
ing interaction in a d-wave superconductor. The local 
pairing interaction was calculated within a real space for- 
mulation of the spin-fluctuation mechanism, and this for- 
malism was used to study the single-impurity problem. 
We find that the local variations in the charge densities 
or hopping integrals may significantly enhance the local 
d-wave gap close to the impurity sites due to a locally en- 
hanced transverse spin susceptibility. We have discussed 
the observation of such enhancement effects in relation 
to STS measurements on the surface of cuprate super- 
conductors. Our model provides a powerful tool to study 
ideas about the effect of inhomogeneity on superconduc- 
tivity at a microscopic level. It will be interesting to find 
whether inhomogeneous defects or magnetic structures 
can enhance the critical temperature, as was found in 
purely phenomenological approaches.— 
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FIG. 6: Real space diagrams for longitudinal spin 
fluctuations to fifth order in U. 



The pairing interaction is dependent only on two spatial 
variables since U is onsite. We get 



V$ (a, c; b, d) = S ab S cd [US ad + U 3 £ Xa 



.Tt 
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+u 5 V y^VtviV 

+ u 2^ X aeX e fXfgXg 
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s ab s cd [ui + u 3 x u x n + u 5 (x u x n f + ---U 

u 3 x u x n 



SabScd Ulad + 



1 - U 2 x u X n 



a d 



(A4) 



Similarly, in the case of the ladder diagrams shown in 
Fig. [7] we get 

V c f (a, c; 6, d) = S ad S bc [U X H + U 3 ^ xltxlt 



Appendix A: Effective pairing interaction in real 
space 

The interaction Hamiltonian is given by 

H- mt = ^s(a,c;b,d)(cl^c di Caic a f + H.c.), (AI) 

a,b,c,d 

where a, ..d are real space points. We define the nonintcr- 
acting spin-dependent susceptibility of the normal phase 
as 



fb 



X a al(T) = -QabAr)Qba^(-T). 



(A2) 



Note the sign convention. The effective interaction in 
real space is evaluated directly from the real space ver- 
sion of the standard spin-fluctuation diagrams. Since the 
Coulomb interaction is local both in space and time, it 
connects fermions at the same site and of opposite spin. 
Omitting the external fermion lines and evaluating the in- 
teraction part from the bubble diagrams shown in Fig. [5] 
gives 



(a, c; b, d) = US ab S cd 5 a d+ 



U 2_, SabbcdGae^Gea^Ged^Gde,^ + 



(A3) 



+u 5 J2xlixlWf 

SadSb c [U\^ + U 3 x n2 + U 5 x^ + 

U 2 x n 
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(A5) 



l-Ux n 

The interaction Hamiltonian now becomes 

H ln t = V oS (°. C ; & ' d ) ( C U C dl C ci c at + H.C.) 

a,b,c,d 

+V$ (a, c; b, d) (4 T 4 ;CH c at + H.c.) 

= v ctt (°> d ' a ' d ) ( c it c rfi c ^ ca t + h - c - 

ad 

+ v $ («. & > «) ( c ! t c ii c &i c at + H - c ) ■ ( A6 ) 

ab 

From this it is apparent that the first part corresponds 
to spin-conserving processes while the second describes a 
spin-flip interaction. 



1. Mean field Hamiltonian 

The result for the effective interaction derived in 
Eq. (|A6[) is treated in a mean-field approach. We consider 
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FIG. 7: Real space diagrams for transverse spin fluctuations 
to third order. 



only the singlet channel. 



- W int + H int 



(A7) 



where we have used that the interactions are symmetric 
with respect to interchange of spatial indices. We define 
two gaps which are respectively symmetric (s) and anti- 
symmetric (a) under interchange of spatial indices 
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eff 



((CjlCif) ± (CjfCil)) , 

{{c 3i c lt ) ± (c jt c 4 }) . (A8) 



For the first term in the interaction Hamiltonian, we get 
by a standard mean field decoupling 



H {1) - 2 V^ ( V r f r,r t 



+ A (i),ij( c i4- C it + c iT c u) + A (i),iMt c U + c 4 c lt) 



The second term gives 



V- 



(2) 



Z) 2 [ A (2).«( c ^t ~ c jT^) + A (2),«(clt c ]4 - c l c it) 



+ A (2),y(c J -;Ci t + c, t c a ) + A^ 2)iij .(cl t ct 4 + c^ct t ) 

Adding these results, we obtain the final effective mean 
field Hamiltonian in the singlet channel 



rrMF 

singlet 2 



where 



AS V eS (i,j) n v 

A ii = 7y {(CjlCit) ~ (Cjt^i))- 



(A9) 



(A10) 



Here, the total effective pairing interaction V e s(i, j) is the 
sum of the effective interactions and V^g derived in 
Eqs. (|A4[) and (|A5|) . Thus, the effective pairing is given 
by the real space matrix 



V e s(i,j) 



v£\i,j) + v£\i,j) 
u 3 x u x n 



Ul 



1 - U 2 x U X n « i - U\ n 



(All) 



which reduces to the usual expression in the translation- 
ally invariant case. 
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